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ACOUSTO-OPTICAL  VECTOR  MATRIX  PRODUCT  PROCESSOR: 
IMPLEMENTATION  ISSUES 


I .  INTRODUCTION 

An  acousto-optic  (AO)  vector  matrix  multiplier  (VMM)  may  be 
conceived  of  as  a  high  speed  pipelined  processor.  Present  AO 
technology  promises  to  allow  multiplication  of  128-element 
vectors  by  a  128  x  128  element  matrix  at  a  rate  of  105  vector- 
matrix  products  per  second.  The  vector  components,  along  with 
the  matrix  elements,  are  complex  numbers  whose  real  and 
imaginary  parts  are  each  8-bit  signed  numbers.  While  such  a 
processor  may  be  used  for  any  linear  transformation,  our 
immediate  interest  is  in  the  computation  of  discrete  Fourier 
transforms. 

The  rate  of  105  such  vector-matrix  products  per  second 
requires  eight-bit  integer  products  at  a  rate  of  approximately 

I. 7  x  1012  multiplies  per  second.  This  rate  is  not  achievable 
with  current  digital  technology.  A  prototype  processor  (rack 
mountable  19“  chassis  not  to  exceed  9  inches  in  height)  using  AO 
technology  was  designed  and  assembled  by  Westinghouse  R&D  Center 
and  evaluated  by  the  Applied  Optics  Branch  of  the  Naval  Research 
Laboratory.  This  report  consists  of  a  descriptive  system 
overview,  a  component  by  component  performance  evaluation,  a 
system  performance  evaluation,  and  suggestions  for  future 
development  in  this  area. 

II.  SYSTEM  OVERVIEW 

A  schematic  diagram  of  the  VMM  processor  is  illustrated  in 
Figure  1.  The  rf  power  driving  the  transducer  of  the  Bragg  cell 
is  modulated  by  a  digitally  controlled  network  of  attenuators. 

The  real  and  imaginary  parts  of  each  component  of  the  input 
vector  are  stored  in  a  buffer  on  the  digital  control  board.  On 
this  digital  control  board  the  data  are  modified  to  correct  for 
the  acoustic  attenuation  in  the  Bragg  cell  and  for  non- 
uniformities  in  the  transverse,  horizontal  intensity  profile  of 
the  laser  beam.  The  data  are  then  clocked  out  to  the  attenuator 
network  at  25.6  MHz.  Immediately  following  the  downloading  of 
the  last  vector  component,  the  information  for  the  entire  vector 
is  represented  by  the  longitudinal  distribution  of  acoustic  power 
in  the  Bragg  cell.  The  laser  is  fired,  and  the  vector  becomes 
represented  by  the  transverse,  horizontal  intensity  distribution 
of  the  diffracted  radiation.  This  diffracted  radiation, 
following  collimation,  is  incident  on  a  rectangular  array  of 
apertures.  The  area  of  each  aperture  represents  the  real  or 
imaginary  part  of  a  matrix  element.  The  radiation  transmitted 
through  the  aperture  array  mask  is  directed  onto  a  linear 
detector  array  such  that  the  radiation  from  each  row  of  apertures 
falls  on  a  single  pixel  of  the  detector  array.  The  output  vector 
is  read  out  of  the  detector  array,  digitized,  and  uploaded  to  the 
host  computer.  In  the  host  computer  the  data  may  be  processed  to 
correct  for  noise  associated  with:  pixel-to  pixel-variations  in 
detector  sensitivity,  detector  dark  current,  and  for  non- 
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uniformities  in  the  vertical  intensity  distribution  of  the 
radiation  field. 

III.  PERFORMANCE  ANALYSIS 

A.  Laser  Performance 

The  laser  is  a  Multiheterostructure  GaAlAs  laser  diode 
operating  at  890  nm  with  a  minimum  rated  peak  radiant  power  of  12 
Watts  (at  300  K)  .  The  exit  pupil  of  the  laser  diode  is  400  fim  by 
2  /nm.  The  laser  is  mounted  so  that  its  longitudinal  axis  is 
horizontal  and  the  long  axis  of  the  exit  pupil  is  vertical.  This 
assures  that  the  laser  beam  propagates  horizontally  with  the 
greatest  beam  divergence  in  the  horizontal  plane.  The  system 
requires  a  laser  that  meets  the  following  performance  criteria: 

energy  per  pulse  -sufficient  to  saturate  the  photodiode  array 

after  passing  through  the  system. 

duration  of  pulse  -sufficiently  short  to  freeze  the  motion  of 

the  acoustic  field  in  the  Bragg  cell  during  a 
laser  firing. 

repeatability  -there  are  no  provisions  to  correct  for  pulse 

to  pulse  fluctuations  in  the  beam  profile. 
Thus,  both  the  horizontal  and  vertical  beam 
intensity  profile  must  be  highly  repeatable. 

firing  rate  -fast  enough  to  permit  the  laser  to  fire  every 

10  /ns,  the  time  needed  to  load  a  vector  into 
the  Bragg  cell. 

The  relative  horizontal  intensity  profile  of  the  laser  was 
measured  using  a  1024-element  photodiode  array  (Reticon  RL1024S) . 
A  typical  horizontal  intensity  profile  is  shown  in  Fig.  2.  The 
VMM  uses  the  portion  of  this  profile  that  extends  approximately 
112  pixels  to  either  side  of  the  beam  axis.  At  the  edges  of 
this  region  the  intensity  dropped  to  88  percent  of  the  peak 
(center)  value.  At  any  pixel  of  the  array,  the  pulse-to-pulse 
variation  of  the  intensity  was  measured  and  found  to  be  approxi¬ 
mately  1%. 

The  full  width  at  half  maximum  (FWHM)  of  the  laser  pulse 
(Fig.  3)  is  10  ns,  less  than  half  of  the  39  ns  (164  /n m )  feature 
size  in  the  Bragg  cell  (see  Section  III.  B.).  Since  one-to-one 
imaging  of  the  Bragg  cell  onto  the  aperture  array  mask  was  used, 
the  horizontal,  center-to-center  aperture  spacings  were  made 
equal  to  the  feature  size  (164  /nm)  in  the  Bragg  cell  (Fig.  4). 
During  the  20  ns  (see  Fig.  3)  duration  of  the  laser  pulse  the 
features  of  the  acoustic  field  in  the  Bragg  cell,  and  thus  the 
features  of  the  diffracted  radiation,  travel  about  82  /nm,  half 
the  center-to-center  aperture  spacing.  To  prevent  effects  of 
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this  motion  from  reaching  the  detector,  the  maximum  width  of  an 
aperture  is  limited  to  82  /Ltm. 

A  discussion  of  the  power  requirements  placed  on  the  laser 
in  connection  with  the  detector  array  performance  can  be  found  in 
Section  III.  D. 

B.  Bragg  Cell,  RF  Attenuator  Network,  and  Digital  RF  Controller 

The  information  is  encoded  in  the  radiation  field  via  an 
Interaction  Model  AOD150  dense  flint  glass  Bragg  Cell.  The  cell's 
specifications  are  listed  in  the  table  below. 


BRAGG  CELL  SPECIFICATIONS 


entrance  pupil 
center  frequency 
acoustic  speed 
measured  attenuation 
(at  150  MHz) 
bandwidth 

time  bandwidth  product 


45  mm  x  3  mm 
150  MHz 
4100  m/s 
«2 . 3  dB/cm 

100  MHz 
1000 


The  power  level  of  the  150  MHz  rf  signal  that  drives  the 
transducer  is  determined  by  the  state  of  the  digitally  controlled 
rf  attenuator  network  (see  Fig.  5) .  Since  the  output  of  the 
digital  controller  is  clocked  at  25.6  MHz,  in  10  ms  the  rf  power 
level  is  specified  256  times  (at  39  ns  intervals) .  Since  the 
speed  of  sound  in  dense  flint  glass  is  4100  m/s,  the  acoustic 
field  generated  in  a  10  /is  interval  is  distributed  over  a  4.1  cm 
length  of  the  cell.  The  digital  controller  thus  partitions  the 
acoustic  field  into  256  sections  (features),  each  section  being 
164  Mm  long.  The  instantaneous  acoustic  power  in  any  particular 
feature  is  the  product  of  the  acoustic  power  output  of  the 
transducer  when  that  feature  was  clocked,  and  the  attenuation  of 
the  signal  after  propagating  from  the  transducer  to  its  present 
position  of  the  feature.  The  intensity  of  the  optical  radiation 
diffracted  by  the  acoustic  field  varies  with  acoustic  power. 

Thus  the  horizontal  intensity  profile  of  the  diffracted  optical 
field  is  also  partitioned  into  256  sections. 

The  150  MHz  rf  signal  driving  the  Bragg  cell  transducer  is 
gated  through  the  attenuator  network  (see  Fig.  5)  by  eight 
digitally  controlled  switches.  The  status  (on/off)  of  each 
switch  is  determined  by  a  digital  control  line.  The  controller 
receives  data  and  control  information  from  an  IBM  PC-AT  via  a 
National  Instruments  GPIB  interface.  The  data  consist  of  two 
vectors  (arrays) ,  VectorData  and  Correct ionData .  These  arrays 
are  stored  in  high  speed  memory  buffers  which  reside  on  the 
digital  control  board.  VectorData  is  the  VMM  input  vector.  In 
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order  to  correct  for  both  the  acoustic  attenuation  within  the 
Bragg  cell  and  nonuni formities  in  the  horizontal  intensity 
profile  of  the  incident  radiation,  the  elements  of  VectorData  are 
multiplied  by  the  corresponding  elements  of  CorrectionData .  For 
this  purpose  the  controller  has  a  resident,  high  speed  8x8  bit 
multiplier  (MC10901) .  Each  byte  of  the  output  of  the  digital 
control  board  is  the  high  order  byte  of  the  product  of 
corresponding  elements  of  VectorData  and  CorrectionData.  Each 
bit  of  these  output  bytes  determines  the  state  of  the 
corresponding  rf  switch  in  the  attenuator  network. 

The  diffraction  efficiency  of  the  Bragg  cell  varies  as  the 
square  of  the  sine  of  the  square  root  of  the  acoustic  power1. 

For  this  particular  Bragg  cell  driven  at  159  MHz  this  relation  is 


where  1^  is  the  diffracted  intensity,  I,-,  is  the  incident 
intensity,  and  P  is  the  acoustic  power  m  watts.  The  prototype 
VMM  is  restricted  to  operate  in  the  domain  where  the  diffraction 
efficiency  varies  linearly  with  acoustic  intensity.  If  eight  bit 
accuracy  is  to  be  maintained,  this  restriction  limits  the 
acoustic  power  to  16  mW  and  the  diffraction  efficiency  to  1.2%. 
Higher  diffraction  efficiencies  might  be  achieved  by  driving  the 
Bragg  cell  at  higher  power  levels  while  correcting  for  the 
nonlinearities.  These  corrections  would  involve  preprocessing 
the  signal  using  a  high  speed  look  up  table.  In  the  VMM 
prototype  this  enhancement  was  not  implemented. 

The  attenuation  network  consists  of  eight  rf  attenuators 
arranged  in  parallel  branches  (Fig.  6) .  The  attenuation  in 
branch  i+1  is  6  dB  greater  than  that  in  branch  i,  i  =  1,..,7. 
Consequently,  the  output  voltage  of  branch  i+1  is  half  that  of 
branch  i.  When  the  voltages  are  recombined  (in  phase),  the 
voltage  of  the  combined  signal  is  the  sum  of  the  voltages  of  the 
branch  signals.  Thus  the  voltage  of  the  network  output  signal  is 
K'V^,  where  the  integer  K  =  0,...,255,  and  V^  is  the  output 
voltage  of  the  branch  with  the  greatest  attenuation.  And  thus 
the  power  level  of  the  network  output  signal  is  K2P^  where  P]^  is 
the  power  output  of  the  branch  with  output  voltage  V1# 

As  indicated  above,  it  is  the  voltages  (amplitudes)  that  add 
when  the  branch  signals  combine.  Since  power  is  proportional  to 
the  square  of  the  voltage,  equal  increments  of  voltage  imply 
unequal  increments  in  power.  Ideally,  it  is  the  increments  in  r:: 
power  that  should  be  equal.  However,  it  does  not  appear 
feasible  to  configure  eight  branches  of  attenuators  to  produce 
256  equally  spaced  power  levels. 

In  the  Bragg  cell,  the  acoustic  power  is  attenuated  at 
2.3  dB  per  cm.  Thus  the  attenuation  over  the  4.1  cm  active 
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length  of  the  cell  is  9.7  dB  (90%).  When  the  elements  of  the 
array  VectorData  have  identical  values,  the  horizontal  intensity 
profile  of  the  diffracted  radiation  field  should  be  uniform.  For 
a  uniform  input  vector,  the  elements  of  CorrectionData  must  cause 
a  monotonic  increase  in  the  attenuation  of  the  network,  with  a 
maximum  attenuation  of  9.7  dB.  To  accomplish  this  the  elements 
of  CorrectionData  must  vary  from  255  to  86  according  to  the 
relation 


K[I]  =  255 • 10-0 • 05$VI/f  . 

In  this  relation  the  attenuation  coefficient  (3  is  2.3  dB/cm,  the 
sound  velocity  V  is  4100  m/s,  the  clock  counter  I  is  the  output 
of  a  MOD  256  counter,  and  the  clock  frequency  f  is  25.6  MHz. 

Since  the  K[I]'s  are  integers,  the  values  obtained  from  the 
expression  must  be  appropriately  rounded.  The  value  of  K[255] 
is  86. 

Consider  the  case  when  VectorData [I]  =  255  (I  =  0,...,255) 
and  the  elements  of  CorrectionData  are  optimally  selected  to 
produce  an  acoustic  field,  and  thus  a  diffracted  infrared  beam, 
with  a  uniform  intensity  profile.  The  distribution  of  acoustic 
power  at  the  time  the  laser  is  fired  can  be  found  in  column  tnree 
of  Table  I.  This  distribution  has  a  mean  of  7476  and  a  standard 
deviation  of  41.  Thus  the  7476  units  of  available  acoustic  power 
can  be  divided  into,  at  most,  91  (7476/82)  distinct  values,  each 
value  an  integral  multiple  of  82  units.  However,  all  91  of  these 
values  cannot  be  realized  because  the  attenuator  network  does  not 
have  sufficient  resolution.  The  attenuator  network  outputs  7396 
(862)  units  of  power  when  the  datum  written  to  it  is  86  and  7225 
(852)  units  when  the  datum  is  85.  Thus  the  minimum  change  in 
power  that  the  attenuator  network  can  deliver  is  171  (7396-7225) 
units.  This  implies  that  the  7476  units  of  available  acoustic 
power  can  be  divided  into  no  more  than  43  (7476/171)  values. 
Consequently,  there  is  a  corresponding  resolution  to  the 
intensity  level  of  the  diffracted  optical  radiation. 

Let  1^  denote  the  optical  intensity  diffracted  by  the 
segment  of  the  Bragg  cell  nearest  the  transducer  when  the  rf 
power  driving  the  transducer  is  The  diffracted  optical 

intensity  level  Id  may  range  from  zero  to  862Ii  while  the 
resolution  of  this  intensity  level  is  1711^.  To  provide 
significantly  distinct  levels  of  diffracted  intensity,  the  rf 
power  must  be  incremented  by  integral  multiples  of  17 IP!.  Thus 
the  diffracted  intensity  level  must  lie  between  zero  and  862Ilf 
and  must  be  an  integral  multiple  of  1711^.  That  is: 

Id  =  N* 171IJ  ±  861!,  N=0,...,43. 

With  only  44  distinct  levels  of  diffracted  intensity,  the 
elements  of  the  VMM's  output  vector  are  limited  to  less  than  6- 
bit  accuracy.  In  order  to  achieve  the  specified  8-bit  accuracy, 
a  full  256  significantly  distinct  levels  of  diffracted  intensity 
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must  be  available.  Increases  in  the  number  of  distinct  levels 
of  diffracted  intensity  can  be  achieved  by:  substituting  a  Bragg 
cell  with  less  attenuation,  controlling  the  rf  power  in  a  manner 
that  results  in  rf  power  levels  with  equal  increments,  increasing 
the  number  of  available  rf  power  levels,  and/or  fabricating  an 
infrared  filter  with  a  spatial  variation  in  transmitted  intensity 
appropriate  to  correct  for  the  effects  of  acoustic  attenuation 
within  the  Bragg  cell. 

The  horizontal  intensity  profile  of  the  radiation  incident 
on  the  Bragg  cell  varies  by  about  12%  over  the  4.1  cm  length  of 
the  cell.  This  variation  is  corrected  for  on  the  rf  control 
board  by  appropriately  modifying  the  elements  of  CorrectionData. 
The  additional  rf  attenuations  used  to  provide  this  correction 
are  relatively  small  and  will  not  require  CorrectionData  elements 
less  than  86. 

It  should  be  noted  that  the  attenuation  board  outputs  256 
equal  increments  of  voltage  rather  than  power.  In  this  regime, 
the  power  for  the  maximum  output  voltage  is  2 16  times  larger 
than  the  smallest  increment  between  power  levels  (which  is  equal 
to  the  smallest  nonzero  power  level) .  Thus  the  attenuation  board 
must  be  calibrated  to  an  accuracy  of  one  part  in  216. 

C.  Aperture  Array  Mask 

The  product  of  two  complex  numbers  X  and  A  can  be  imple¬ 
mented  via  matrix-vector  products  involving  only  real  numbers. 


Let 

A  =  aR  +  ia1, 

(1) 

and 

X  =  xR  +  ix1. 

(2) 

Then 

r=A*x^ 

(aRxR  -  a1^1)  +  i(aRxI  +  a-*-xR)  =  pR  +  ip1, 

(3) 

or, 

In  the  VMM,  all  complex  multiplies  are  executed  with  this 
algorithm.  Each  complex  element  of  the  vector  is  represented  as 
a  two-element  real  vector  and  each  complex  matrix  element  is 
represented  as  a  2  x  2  real  matrix.  Thus  a  128  complex-element 
vector  is  represented  by  a  256  real-element  vector  and  a 
128  x  128  complex-element  matrix  is  represented  as  a  256  x  256 
real-element  matrix. 
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When  a  pulse  of  radiation  of  uniform  intensity  I  is  incident 
on  an  aperture  of  area  S,  the  transmitted  energy  per  pulse  is 
proportional  to  the  product  I*S.  In  the  VMM,  the  matrix  is 
represented  as  a  rectangular  array  of  256  x  (256  +  1)  apertures 
(see  Fig.  7) .  The  area  of  an  aperture  determines  the  value  of 
the  real  number  it  represents.  An  offset  is  used  to  represent 
signed  numbers:  aperture  areas  less  then  a  specified  aperture 
size  SD  represent  negatively  valued  numbers  while  aperture  areas 
greater  than  SQ  represent  positively  valued  numbers.  The  first 
row  of  apertures,  each  of  area  SQ,  allows  calibration  for  the 
bias  or  offset  associated  with  the  choice  of  SQ.  In  the 
postprocessing  of  the  electronic  output,  all  levels  are  compared 
to  this  bias  level.  The  input  vector  is  represented  by  the 
transverse,  horizontal  intensity  distribution  of  the  radiation 
field  incident  on  the  aperture  array.  This  horizontal  intensity 
distribution  is  divided  into  256  segments  (features) ,  each 
segment  incident  upon  a  single  column  of  apertures.  The 
intensity  level  of  a  segment  corresponds  to  the  value  of  the 
number  it  represents.  A  specified  intensity  I0  corresponds  with 
the  value  zero,  while  intensities  less/more  than  IQ  correspond 
to  negatively/positively  valued  numbers. 

A  vector  matrix  product  can  be  represented  as  fellows: 

Yi  =  ?  aij  Xj  (5) 

where  the  quantities  are  S-bit  signed  integers. 

Let  a^  '  =  a^j  +  Bi  ,  Xj  '  =  Xj  +  B2  ,  and 

y i 1  =  ?  (aij  +  ®l)  (Xj  +  B2 )  •  (6) 

To  insure  that  the  values  of  the  primed  parameters  defined  above 
are  non-negative,  the  values  of  the  bias  terms  Bi  and  B2  are  both 
set  to  127.  The  values  of  these  primed  quantities  are 
represented  by  aperture  area  and  radiation  intensity.  Combining 
expressions  (5)  and  (6)  yields  a  relation  expressing  yi  in  terms 
of  yi'  and  a  sum  involving  three  bias-dependent  terms. 

Yi  =  Yi'  -  ?  (ai  j  B2  +  Xj  +  B1B2)  (7) 

In  order  to  determine  the  values  of  the  voltages  associated  with 
B^  and  B2  we  input  the  zero  vector  (x-}  =  0  for  all  j)  in  the  VMM. 
With  a  zero  input  vector  the  intensity  of  the  radialion  incident 
on  each  column  of  the  aperture  array  is  IQ  (corresponding  to  B2) . 
The  radiation  from  the  bias  row  of  apertures  is  collected  at  the 
zeroth  pixel  of  the  detector  array.  The  detector  output  voltage 
Vi  associated  with  the  charge  integrated  by  this  radiation  during 
one  laser  pulse  corresponds  to  the  value  of  the  term  y0  =  0 . 

This  voltage  also  corresponds  to  the  value  of  the  term 

V1  =  Yo '  =  N*BiB2,  where  N  (=256)  is  the  number  of  apertures  in 
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the  bias  row.  This  zero  vector,  when  applied  to  the  bias 
calibration  row  at  the  mask  (i=0)  yields  the  term  _  N*B^B2, 
where  N  (=256)  is  the  number  of  elements  in  a  row.  Since  voltage 
levels  corresponding  with  both  B^  and  B2  are  needed,  a  second 
calibration  measurement  is  required  to  specify  voltage  levels 
corresponding  to  both  3^  and  B2.  This  is  done  by  multiplying  the 
matrix  by  a  vector  whose  elements  are  all  equal  to  the  maximum 
value  possible.  That  is,  Xj  =  128  for  all  j.  Then  Eq.  (7) 
yields  V2  r;  yo'  *  N*B^(128  +  B2).  Simultaneously  solving  for  the 
bias  tenas  results  in  the  relations 

R  =  V2  ~  Vj  _  _  128V]_ 

Bl  128N  V2  -  Vx  • 

As  each  aperture  represents  an  8-bit  signed  number,  there 
are  256  distinct  aperture  sizes  allowed.  The  architecture  is 
such  that  each  aperture  has  a  height  of  50  urn  and  a  width  that 
may  vary  between  zero  and  82  Mm  (Fig.  7)  .  The  optical  power 
passing  through  an  aperture  is  proportional  to  the  product  of  the 
width  of  the  aperture  and  the  intensity  of  the  incident 
radiation.  However,  due  to  diffraction,  not  all  the  transmitted 
power  is  delivered  to  the  detector.  The  aperture  widths  are 
selected  so  that  when  the  transverse  intensity  distribution  of 
the  incident  radiation  is  uniform  and  constant,  the  power 
distribution  delivered  to  the  detector  via  each  aperture  will  be 
one  of  256  values  which  are  equally  spaced. 

A  cylindrical  and  a  spherical  lens  combination  (Fig.  8) 
collects  the  radiation  from  the  aperture  array  and  directs  it  to 
the  linear  detector  array.  With  this  arrangement  the  detector 
array  is  positioned  in  the  image  plane  of  the  aperture  array 
mask  vertically  and  the  Fourier  transform  plane  of  the  aperture 
array  mask  horizontally.  Consider  the  optical  power  delivered  to 
a  detector  pixel  from  a  single  aperture.  Assuming  the  aperture 
to  be  uniformly  illuminated,  the  horizontal  intensity  distri¬ 
bution  (Fig.  9)  at  the  detector  pixel  is  proportional  to  the 
square  of  a  sine  function 

Relative  Intensity (x)  =  sine2 (  x/x')  (8) 

where  x'  =  (lower  case  lambda)  f/w,  lambda  being  the  wavelength, 
f  the  focal  length  and  w  the  aperture  width.  The  incident 
relative  power  delivered  to  the  pixel,  the  integral  of  the 
relative  intensity  over  the  full  2500  Mm  width  of  the  detector 
pixel,  varies  with  aperture  width  w.  The  aperture  widths  wn 
(n=0, . . . , 255)  are  selected  to  produce  a  relative  power  spectrum 
Pn  =  nP0/255,  where  PQ  is  the  relative  power  incident  on  the 
pixel  when  the  aperture  width  w  is  the  maximum  82  Mm.  This 
relation  is 


Pn  =  nP0/255 


P+1250  Mm 

sinc2(  x/xn')  dx 

-1250  Mm 


xn,=f/wn 


(9) 
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The  relative  power  constant  Pq  is  determined  by  integrating 
Eq.  (9)  when  n  is  255  and  wn  is  82  /im.  The  values  of  the 
remaining  wn's  (n=0 , . . . , 254)  are  determined  numerically  from 
Eq.  (9). 

An  aperture  array  mask  based  on  the  above  calculations  and 
designed  for  calculating  discrete  Fourier  transforms  was 
fabricated  using  electron  beam  lithography.  The  output  vector  of 
the  VMM,  when  used  with  this  mask,  is  the  discrete  Fourier 
transform  of  the  input  vector. 

D.  Detector  Array 

The  photodetector  consists  of  a  linear  1024  element 
photodiode  array  (Reticon  RL1024S)  with  each  photodiode  measuring 
25  /jm  by  2500  jum .  The  array  is  configured  to  be  read  out  as  512 
separate  elements,  each  element  consisting  of  an  adjacent  pair  of 
pixels.  The  maximum  readout  rate  of  the  array  is  5  MHz.  The 
widths  and  the  center-to-center  spacing  of  the  pixels  of  the 
photodiode  array  influence  the  design  of  the  aperture  array  mask 
as  well  as  the  supporting  aperture-to-detector  optics.  In 
addition  to  the  bias  row,  the  aperture  array  mask  contains  256 
rows  of  apertures.  The  radiation  from  each  row  of  apertures 
falls  on  a  unique  detector  pixel  pair.  To  prevent  errors  due  to 
crosstalk,  between  each  of  these  '•active"  pixel  pair  is  a  "dummy" 
pixel  pair  whose  output  is  ignored.  The  radiant  energy  per  pulse 
incident  on  an  active  pixel  pair  corresponds  to  an  8-bit  integer, 
the  real  or  imaginary  part  of  an  element  of  the  output  vector. 

The  array-to-detector  optics  were  designed  to  have  a 
magnification  of  unity,  thus  the  aperture  array  masks  are 
designed  so  the  center-to-center  spacing  of  the  rows  is  100  fim, 
twice  the  center-to-center  pixel  pair  spacing  of  the  detector. 

The  saturation  energy  per  unit  area  of  the  Reticon  RL1024S 
is  50  nJ/cm2.  In  order  to  utilize  the  detector's  full  dynamic 
range,  the  maximum  integrated  intensity  per  laser  pulse  should  be 
sufficient  to  drive  it  to  saturation.  The  equation  describing 
this  is: 


Saturation  Energy 
nl  n2  n3  pw 


where  PW  is  the  pulse  width  of  the  laser  and  n^,  r\2 ,  and  n3  are 
the  optical  efficiencies  for  the  Bragg  cell,  the  lenses  and 
mirrors,  and  the  aperture  array  mask.  The  maximum  pulse  width  is 
20  ns  while  n^,  n2 ,  and  n3  are  estimated  to  be  0.012,  0.54,  and 
0.125.  Thus,  to  drive  the  detector  into  saturation  with  a  single 
laser  pulse,  a  minimum  laser  power  of  one  kilowatt  is  required. 
The  supplied  laser  has  a  10  ns  pulse  width  and  12  watts  of  peak 
power  (See  Section  III.  A.),  not  all  of  which  is  incident  upon 
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the  entrance  pupil  of  the  Bragg  cell.  With  the  rf  power  driving 
the  Bragg  cell  at  a  constant  3  W,  we  repeatedly  fired  the  laser 
until  the  detector  became  saturated.  Saturation  was  reached 
after  approximately  100  shots.  The  diffraction  efficiency  at  3  W 
of  rf  power  («1.5  W  of  acoustic  power)  is  about  0.1. 
Unfortunately,  a  maximum  rf  power  of  «32  mW  (16  mW  of  acoustic 
power)  is  necessary  if  the  nonlinear  terms  in  the  relation 
between  the  diffraction  efficiency  and  the  acoustic  power  are  not 
to  detract  from  the  specified  eight  bit  accuracy.  Using  the  12  W 
laser  with  a  pulse  width  of  10  ns,  saturation  should  occur  after 
some  550  laser  firings.  To  saturate  the  detector  with  one  laser 
firing  at  a  pulse  width  of  10  ns  would  require  a  6.6  KW  laser. 

See  the  discussion  at  the  end  of  Section  III.  B.  that  proposes  a 
modification  to  reduce  the  required  laser  power. 

It  takes  10  jus  to  load  a  vector  (256  bytes)  into  the  Bragg 
cell  (and  20  ns  for  the  laser  to  fire) .  Since  this  process  can 
be  repeated  continuously,  a  continuous  serial  transfer  rate  of 
25.6  MHz  seems  feasible.  However,  because  the  maximum  readout 
rate  of  the  detector  is  5  MHz,  the  maximum  continuous  throughput 
of  the  prototype  is  reduced  from  105  vector  matrix  products  each 
second  to  104  products  each  second. 

E.  Imaging  and  Focusing  Optics 

The  optics  that  image  the  mask  onto  the  detector  array 
(Fig.  8)  consist  of  a  cylindrical  lens  and  a  spherical  lens,  each 
with  a  75  mm  diameter  and  a  200  mm  focal  length.  The  cylindrical 
lens  is  positioned,  with  its  axis  horizontal,  one  focal  length 
from  the  mask;  the  spherical  lens  is  positioned  two  focal  lengths 
from  the  mask;  and  the  detector  array  is  located  one  focal  length 
from  the  spherical  lens.  With  this  configuration  the  plane  of 
the  detector  coincides  with  the  image  plane  of  the  mask 
(vertically)  and  with  the  Fourier  transform  plane  (horizontally) . 
Thus  the  radiation  from  each  row  of  the  mask  falls  onto  a  single 
pixel  pair  of  the  detector  array. 

To  perform  acceptably,  the  image  of  each  mask  row  must  fall 
on  the  corresponding  pixel  pair  of  the  detector.  Since  the 
detector  array  contains  512  pixel  pairs,  the  two  lens  combination 
must  have  a  linear  magnification  of  unity  to  one  part  in  512  (or 
0.2%).  This  implies  that  the  focal  lengths  of  the  two  lenses  be 
equal  to  within  0.2%.  In  assembling  the  prototype  no  attempt  was 
made  to  accurately  match  the  focal  lengths.  The  image  of  the 
aperture  array  mask  and  the  pixels  of  the  detector  array  were 
found  to  be  significantly  out  of  registry  (Fig.  10) .  Since  the 
lenses  are  «  f/2.7,  it  is  likely  that  additional  difficulties 
with  registry  can  be  attributed  to  spherical  aberration. 

F.  Data  Acquisition  and  Post  Processing 

In  the  prototype  the  detector  array  is  mounted  on  an  EG&G 
Reticon  (Model  RC1024L)  evaluation  circuit  board.  The  array 
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output  was  digitized  via  a  12-bit  A/D  converter  and  stored  in  a 
buffer  resident  in  the  host  microcomputer.  No  attempt  was  made 
to  process  the  data  to  correct  for  nonuniformities  in  the 
vertical  distribution  of  the  radiation  or  for  pixel-to-pixel 
variations  in  detector  sensitivity.  To  correct  noise  associated 
with  the  response  of  the  detector  to  dark  current,  the  output 
voltage  associated  with  dark  current  was  measured  and  subtracted 
from  the  data. 

IV.  SUMMARY 


In  most  of  our  measurements  the  sequence  of  actions 
implemented  was: 

•  The  input  vector  and  the  correction  vector  were 
downloaded  from  the  microcomputer  to  buffers  on  the  rf 
control  board. 

•  The  detector  array  was  initialized  ( scanned). 1 

•  The  corresponding  elements  of  the  input  vector  and  the 
correction  vector  were  multiplied  and  the  products 
(sequentially)  clocked  into  the  Bragg  cell  via  the 

rf  attenuator  network  and  transducer. 

■  The  laser  was  fired. 

•  The  detector  array  was  read  (scanned) ,  and  the  output 
data  were  digitized  and  stored. 

The  prototype  VMM's  performance  fell  short  of  meeting  its 
specifications.  The  causes  of  the  shortcomings  include: 

•  The  acoustic  attenuation  (2.3  dB/cm)  in  the  dense  flint 
glass  Bragg  cell  is  unacceptably  high. 

•  The  increments  in  rf  power  levels  are  not  equal.  This 
reduces  the  dynamic  range  by  at  least  a  factor  of  two. 

■  Remaining  within  the  linear  region  of  the  diffraction 
efficiency/acoustic  power  curve  places  an  unacceptably  lev 
(1.2%)  upper  bound  on  the  diffraction  efficiency. 

■  The  laser  energy  per  pulse  is  two  orders  of  magnitude  too 
low  to  provide  sufficient  dynamic  range  at  the  detector. 

•  Due  to  poor  optical  design  and  component  quality,  the 
radiation  field  features  were  out  of  registration  with  the 
aperture  and  detector  arrays. 

•  The  detector  readout  rate  (>100Ms/vector)  and  digitization 
rate  are  unacceptably  low,  and  there  are  no  corrections  for 
noise  associated  with:  dark  current,  nonuniformities  in  the 
vertical  distribution  of  the  radiation,  and  variations  in 
sensitivity  from  pixel-to-pixel. 

Suggested  remedies  for  some  of  these  shortcomings  are: 

•  Replace  the  dense  flint  Bragg  cell  with  one  made  of  Te02 . 

The  substitution  of  a  longitudinal  Te02  Bragg  cell  for  the 
dense  flint  glass  cell  would  decrease  the  acoustic 


1.  The  array  was  initialized  prior  to  firing  the 
laser  in  order  to  minimize  noise  associated  with  both 
dark  current  and  ambient  radiation. 
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attenuation  by  8  dB  while  increasing  the  diffraction 
efficiency  per  watt  of  acoustic  power  by  a  factor  of  3.8. 

The  acoustic  velocity  in  longitudinal  TeC>2  is  4200  m/s, 
almost  the  same  as  the  4100  m/s  acoustic  velocity  in  dense 
flint  glass. 

•  Generate  the  150  MHz  signal  that  drives  the  Bragg  Cell  with 
a  300+  MHz  12-bit  digital  to  analog  (D/A)  converter.  To 
implement  such  an  rf  source,  a  high  speed  processor  is 
needed  to  process  the  input  vector  data  and  thus  correct  for 
the  nonlinearities  in  the  diffraction  efficiency/acoustic 
power  relation.  A  D/A  converter  meeting  the  above 
specifications  is  beyond  the  reach  of  current  technology. 

•  Select  lenses  with  more  closely  matched  focal  lengths  and 
less  spherical  aberration.  Adequate  corrections  for 
spherical  aberration  will  likely  require  higher  f  numbers 
(and  longer  focal  lengths) . 

■  Fabricate  a  detector  array  with  a  faster  readout  rate.  Use 
real  time  hardware  rather  than  post  processing  software  to 
correct  for  noise  associated  with:  dark  current, 
nonuniformities  in  the  vertical  distribution  of  the 
radiation,  and  pixel-to-pixel  variations  in  sensitivity. 

■  Obtain  a  laser  source  with  more  power  and  less  pulse-to- 
pulse  variation  in  power  and/or  beam  profile.  If  laser 
diode  arrays  continue  to  improve  they  may  prove  to  be 
suitable. 
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TABLE  I 


CORRECTION  FOR  THE  ATTENUATION  OF  ACOUSTIC  POWER 


CLOCK 

INJECTED 

ATTENUATED 

DEVIATION 

COUNT 

POWER* 

POWER* 

I 

K2[I] 

PA** 

(PA2-<PA2> 

0 

65025 

7478 

2 

1 

64516 

7482 

6 

2 

64009 

7487 

11 

3 

63504 

7491 

15 

4 

63001 

7495 

19 

5 

62500 

7499 

23 

6 

62001 

7502 

26 

7 

61504 

7505 

29 

8 

60516 

7448 

-28 

9 

60025 

7450 

-26 

10 

• 

59536 

• 

7452 

-24 

• 

• 

247 

• 

7921 

7401 

• 

-75 

248 

7921 

7464 

-12 

249 

7921 

7528 

52 

250 

7744 

7422 

-54 

251 

7744 

7486 

10 

252 

7744 

7549 

73 

253 

7569 

7442 

-34 

254 

7569 

7505 

29 

255 

7396 

7396 

-80 

MEAN  7476 
STANDARD  DEVIATION  41 

* 

** 


Numbers  represent  relative  values. 
PA  =  K2[I] •10"°'1^v/f 
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I _ I 


OPTICAL  SYSTEM 


Fig.  1  —  Schematic  of  system  including  the  host  computer,  laser,  Bragg  cell,  aperture  array  mask, 
detector,  rf  power  supply,  rf  attenuator  network,  digital  control  board,  and  the  digital  I/O  &  A/D 
converter  board. 


Fig.  2  —  The  horizontal  intensity  profile  of  laser  beam  measured  a  distance  of  45  cm  in  front  of  the 
laser.  Only  the  shaded  region  falls  within  the  entrance  pupil  of  the  VMM’s  optical  system. 


Fig.  3  —  The  intensity  of  the  laser  pulse  versus  time.  The  FWHM 
of  the  pulse  is  ~  10  ns. 
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DIGITAL  CONTROL  BOARD 


Fig.  5  —  Schematic  of  the  Bragg  cell,  the  rf  attenuator  network,  and 
the  digital  control  board. 


Fig.  6  —  The  rf  attenuator  network  consists  of  eight  branches,  each  successive  branch  having  6  dB 
additional  attenuation.  Each  branch  also  contains  a  two  state  (on/off)  digitally  controlled  switch  and 
an  adjustable  phase  shifter  (not  shown). 


i  50  (im  high  with  widths  up  to  82  ^m.  The  center-to-center  spacings  are  100  vertically  and 
fim  horizontally. 


OPTICAL  LAYOUT 
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Fig.  8  —  Overview  of  the  optical  layout.  The  spherical  lenses  are  denoted  with  L’s  and  the  cylindri¬ 
cal  lenses  with  C’s.  The  lenses  all  have  focal  lengths  of  200  mm  which  is  the  distance  between  the 
vertical  lines  on  the  grid  to  the  right  of  C,.  There  are  front  surface  mirrors  at  the  dotted  lines,  each 
deflecting  the  beam  by  90°. 


Fig.  10  —  The  readout  of  the  array,  displayed  on  an  oscilloscope  screen,  shows  that  the  registration 
of  aperture  row  images  and  every  other  detector  pair  was  inadequate.  If  the  registry  was  adequate, 
the  two  curves  would  never  cross. 
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